Quantum Interference in Three Photon Down Conversion 
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We study degenerate three photon down conversion as a potential scheme for generating nonclas- 
sical states of light which exhibit clear signatures of phase space interference. The Wigner function 
representing these states contains an interference pattern manifesting quantum coherence between 
distinct phase space components, and has substantial areas of negativity. We develop an analytical 
description of the interference pattern, which demonstrates how the oscillations of the Wigner func- 
tion are built up by the superposition principle. We analyze the impact of dissipation and pump 
fluctuations on the visibility of the interference pattern; the results suggest that some signatures of 
quantum coherence can be observed in the presence of moderate amount of noise. 

PACS number(s): 42.50.Dv, 03.56.Bz 



I. INTRODUCTION 



The superposition principle is a fundamental ingredi- 
ent of quantum theory, resulting in interference phenom- 
ena not existing in classical mechanics. In atomic, molec- 
ular, and optical physics this striking feature of quantum 
mechanics can be studied within several examples of sim- 
ple quantum systems: a trapped ion, a diatomic molecule 
and a single electromagnetic field mode in a cavity or 
free space. In this context Schleich and Wheeler |]| de- 
veloped a phase space picture of quantum interference. 
They demonstrated in the semiclassical limit, that quan- 
tum mechanical transition probabilities are governed by 
a set of simple rules in the phase space: a probability 
amplitude is given by a sum of overlap areas, with phase 
factors defined by an area caught between the states. 

Recent developments in quantum optics have gener- 
ated significantly increased interest in the phase space 
representation of quantum states, providing feasible 
schemes for measuring the Wigner functions of a single 
light mode ^-0], the vibrational manifold of a diatomic 
molecule ||, and the motional state of a trapped atom 
H , or an atomic beam Q] . These advances open up new 
possibilities in experimental studies of the quantum su- 
perposition principle, as the Wigner function provides 
direct insight into the interference phenomena through 
its fringes and negativities, and also completely charac- 
terizes the quantum state. Additionally, negativity of the 
Wigner function is a strong evidence for the distinctness 
of quantum mechanics from classical statistical theory. 
Consequently, it is now possible to obtain full information 
on the coherence properties of a quantum state by mea- 
suring its Wigner function, instead of observing quantum 
interference only as fringes in marginal distributions of 
single observables. 

Therefore, schemes for generating quantum states 
with nontrivial phase space quasidistributions, especially 
those possessing substantial negativities, are of consider- 
able interest. The system that appears to provide the 
most opportunities currently is a trapped ion, whose 



quantum state can be quite easily manipulated through 
interaction with suitably applied laser beams S. In the 
case of travelling optical fields, the range of available in- 
teractions is far more restricted, and generating states 
with interesting phase space properties is a nontrivial 
task from both theoretical and experimental points of 
view. One of the states that most clearly illustrate quan- 
tum interference is a superposition of two distinct coher- 
ent states whose generation in microwave frequency 
range has been recently reported flo|| . Production of 
these states in the optical domain has been a subject 
of considerable theoretical interest. Though several inge- 
nious schemes have been proposed [ll|-Q, they require 
extremely precise control over the dynamics of the sys- 
tem, which makes them very difficult to implement ex- 
perimentally. 

In this paper we study degenerate three photon down 
conversion |l5|-|2^] as a scheme for generating states of 
light that exhibit clear signatures of phase space inter- 
ference. This generation scheme seems to be quite at- 
tractive, since, as we will show, it is not overly sensitive 
to some sources of noise. Additionally, numerous exper- 
imental realizations of two photon down conversion for 
generating squeezed light give a solid basis for studying 
higher order processes, at least in principle, and devel- 
opments in nonlinear optical materials suggest it may be 
possible to re-examine higher-order nonlinear quantum 
effects. 

Interference features of states generated in higher or- 
der down conversion have been first noted by Braunstein 
and Caves fl7| l, who showed oscillations in quadrature 
distributions and explained them as a result of coher- 
ent overlap of two arms displayed by the Q function. 
The purpose of the present paper is to provide a detailed 
analysis of the interference features, based on the Wigner 
function rather than distributions of single observables. 
Compared to the Q function, discussed previously by 
several authors, the Wigner function carries explicit ev- 
idence of quantum coherence in the form of oscillations 
and negative areas. These features are not visible in the 
Q function, which describes inevitably noisy simultane- 
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ous measurement of the position and the momentum. 

The states generated in three photon down conversion 
cannot be described using simple analytical formulae. It 
is thus necessary to resort to numerical means in order 
to discuss their phase space properties. However we will 
show that the interference features can be understood 
with the help of simple analytical calculations. These 
calculations will reveal the essential role of the superpo- 
sition principle in creating the interference pattern in the 
phase space. Experimental realization of the discussed 
scheme along with detection of the Wigner function of 
the generated field would be an explicit optical demon- 
stration of totally nonclassical quantum interference in 
the phase space. 

This paper is organized as follows. First, in Sec. |lj], we 
discuss some general properties of the Wigner function. 
In Sec. HI we present the numerical approach used to deal 
with three photon down conversion. The Wigner function 
representing states generated in this process is studied in 
detail in Sec. IV. In Sec. [v| we discuss briefly prospects 
of experimental demonstration of quantum interference 



using the studied scheme. Finally Sec. VI concludes the 
results. 



II. GENERAL CONSIDERATIONS 

Before we present the phase space picture of three pho- 
ton down conversion, let us first discuss in general how 
the interference pattern is built up in the phase space by 
the superposition principle. Our initial considerations 
will closely follow previous discussions of the semiclassi- 
cal limit of the Wigner function [^3| . They will give us 
later a better understanding of the interference we are 
concerned with in three photon down conversion, and 
help us to derive an analytical description of the inter- 
ference pattern for this specific case. 

We will start by considering a wave function of the 
form 



ip(q) = A(q)exp[iS(q)], 



(1) 



where S(q) is a real function defining the phase and A(q) 
is a slowly varying positive envelope. The Wigner func- 
tion of this state is given by (throughout this paper we 
put h — 1): 

W^(q,p) = ^- J dxA(q-x/2)A{q + x/2) 

x exp[— ipx — iS(q — x/2) + iS(q + x/2)]. 

(2) 

Let us first separate the contribution from the direct 
neighborhood of the point q. For this purpose we will 
expand the phase S(q) up to the linear term and take 
the value of the envelope at the point q, which gives 



Thus, this contribution is localized around the momen- 
tum S'(q), and creates a concentration along the "tra- 
jectory" (q,p — S'(q)). This result has a straightforward 
interpretation in the WKB approximation of the energy 
eigenfunctions, where the phase S(q) is the classical ac- 
tion and its spatial derivative yields the momentum. In 
this case, Eq. (||) simply states that the Wigner func- 
tion contains a positive component localized along the 
classical trajectory p3[ . 

We will now study more carefully the relation between 
the wave function and the Wigner function, taking into 
account contributions from other parts of the wave func- 
tion, denoted symbolically by dots in Eq. (||). To make 
the discussion more general, we will take the wave func- 
tion to be a superposition of finite number of components 
defined in Eq. (|l|): 



4>{q) = ^2Ai(q) exp[iSi(q)]. 



(4) 



The Wigner function is in this case a sum of integrals 

w^p) = yYj [ ^Mq-x/2)A.j{q + x/2) 

x exp[— ipx — iSi(q — x/2) + iSj(q + x/2)]. 

(5) 

We will evaluate these integrals with the help of the sta- 
tionary phase approximation. The condition for the sta- 
tionary points is given by the equation 



S' i (q-x/2)+S'Aq + x/2) = 2p, 



(6) 



which has a very simple geometrical intepretation. It 
shows that the contribution to the Wigner function at 
the point (q,p) comes from the points of the "trajecto- 
ries" (qi,pi = S-iq)) and (qj,Pj = S' } (q)) satisfying 



(ft + ?j)/2 = q 
(Pi+Pi)/2 =P, 



(7) 



W^q,p)^A 2 (q)S(p-S'(q)) + 



(3) 



i.e. (<?,p) is a midpoint of the line connecting the points 
(qi,Pi) and (qj,pj). These points may lie either on the 
same trajectory, i.e. i — j or on a pair of different ones. 
In particular, for i — j we get that qi = qj = q is always 
a stationary point for p — S^q), which justifies the ap- 
proximation applied in deriving Eq. (||). At these points 
the second derivative of the phase disappears. There- 
fore we will calculate them separately, using the previous 
method. For the remaining pairs, we expand the phases 
up to quadratic terms and perform the resulting Gaus- 
sian integrals. As before, we neglect variation of the en- 
velopes, taking their values at the stationary points. This 
yields an approximate form of the Wigner function: 
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wv(«,p)«2A ? (?)*(p-^(?)) 

i 



v /^(5. i "(g J )-5j'(g J ))/2 
x exp[ip(g 4 - g,) - iS l {q l ) + iSj(qj)}, (8) 



n.ij 

S'(9i) + S'(9,) = 2p 



where the second double sum excludes the case i — j and 
Qi = Qj = q- 

Thus, the Wigner function of the state defined in 
Eq. (Q) exhibits two main features. The first one is 
presence of positive humps localized along "trajectories" 
(qi,Pi = S'(qi)). Any pair of points on these trajectories 
gives rise to the interference pattern of the Wigner func- 
tion at the midpoint of the line connecting this pair. Let 
us note that the result that the interference pattern in 
a given area is generated by equidistant opposite pieces 
of the quasidistribution corresponds to the phase space 
picture of the superposition of two coherent states , 
for which the interference structure lies precisely in the 
center between the interfering states. 



III. NUMERICAL CALCULATIONS 

Numerical results presented in the following parts of 
the paper are obtained using a model of two quantized 
light modes: the signal and the pump, coupled by the 
interaction Hamiltonian: 



H = i\[b{a 



^3 



(9) 



where A is the coupling constant, and a and b are the 
annihilation operators of the signal and pump mode, re- 
spectively. This Hamiltonian is very convenient for nu- 
merical calculations, as it commutes with the operator 
N = 3a^a + tfb, and can be diagonalized separately in 
each of the finite-dimensional eigenspaces of N. Details 
of the basic numerical approach to these kinds of Hamil- 
tonians can be found for example in Refs. JlJ,[li|. In 
contrast, the limit of a classical, undepleted pump is dif- 
ficult to implement numerically due to singularities of the 
evolution operator on the imaginary time axis p5[ . 

We assume that initially the signal mode is in the vac- 
uum state |0), and the pump is in a coherent state \fj). 
After evolution of the system for the time t, which we cal- 
culate in the interaction picture, we obtain the reduced 
density matrix of the signal field by performing the trace 
over the pump mode: 



p(t) = Tr pump [ 



iHt 



|0)(0|® |/3>(/J|e 



(10) 



In general, p(t) describes a mixed state, as the interacting 
modes get entangled in the course of evolution. This den- 
sity matrix is then used to calculate the Wigner function 
and other obervables of the signal mode studied further 
in the paper. In the discussions, we will make use of the 



analogy between a single light mode and a harmonic os- 
cillator, assigning the names position and momentum to 



the quadratures q = (a 
respectively. 



)/v / 2 and p={a- a?)/V2i, 



IV. WIGNER FUNCTION 

We will restrict our studies to the regime of a strong 
pump and a short interaction time. This regime is the 
most reasonable one from experimental point of view, 
as strong pumping allows us to compensate for the usu- 
ally weak effect of nonlinearity, and the short interaction 
time gives us a chance to ignore or to suppress dissipa- 
tion. We can gain some intuition about the dynamics of 
the system by considering the classical case; this is done 
in Appendix [X| The most important conclusion is that 
in the classical picture the origin of the phase space is an 
unstable fixed point, with three symmetric directions of 
growth, in a star-like formation. 

In Fig. |l| we depict the Wigner function representing 
the state of the signal field generated for the parameters 
(5 = 10 and t — 0.025/ A. This state is almost pure, as 
Tr[p 2 ] = 0.92 indicates little entanglement between the 
pump and the down-converted mode. The three devel- 
oping arms follow the classical directions of growth from 
the unstable origin of the phase space. The coherence 
between these components results in an interference pat- 
tern filling the regions between the arms, consisting of 
positive and negative strips. Thus, the Wigner function 
is "forced" by the superposition principle to take nega- 
tive values in order to manifest the quantum coherence 
of the state. 

Let us now study in more detail how the interference 
pattern is generated by coherent superposition of distinct 
phase space components. We will focus our attention on 
the three arms displayed by the quasidistribution, ne- 
glecting the bulk of positive probability at the origin 
of the phase space remaining from the initial vacuum 
"source" state. As the purity factor of the generated 
state is close to one, we will base our calculations on 
pure states. The relation between the wave function and 
the Wigner function derived in Sec. [n] suggests that the 
arms can be modelled by three components of the wave 
function: 



<p(q) = <po(q) + <pi(q) + f2{q) 



Hi) 



with slowly varying envelopes and the position- 
dependent phase factors: So(q) = 0, 5i(g) — \/~3q 2 /2, 
and 52(g) = — \/3g 2 /2, respectively. The interference 
pattern observed in the phase space is a result of the 
coherent superposition of these three components. 

However, in order to calculate quantitatively the struc- 
ture of the interference pattern, we need to know the 
relative phase factors between the wave functions in 
Eq. (pi]). We will obtain these factors with the help of 
the additional information that the Hamiltonian defined 
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in Eq. (J9|) excites or annihilates triplets of signal photons. 
Consequently, the photon distribution of the generated 
state is nonzero only for Fock states being multiples of 
three, as the initial state was the vacuum. Using this fact, 
we can define an operator which performs a rotation in 
phase space by an angle 6: 

U(0) = exp(-i0a f a) (12) 

and impose the relations ipi — U(2Tr/3)ipo and if2 = 
C r (47r/3)< y 9o- This choice for the phase of the operator 
U (9) ensures that the superposition defined in Eq. (|ll| ) 
has the necessary property to generate the correct triplet 
photon statistics. Let us now assume that ipo is given 
by a slowly varying positive function A{q), localized for 
q > 0. We will not consider any specific form of the 
envelope A{q), as the main purpose of this model is to 
predict the position and shape of the interference fringes. 
The other two wave functions can be calculated with the 
help of the formula derived in Appendix [b], which finally 
yields: 

(fo{q) = A(q), 



W int (q,p) = 3l/4?rl/2 



As the envelope A{q) is a positive function, the oscilla- 
tions of the interference pattern are determined by the 
argument of the cosine function. The lines of constant 
argument are hyperbolas with asymptotics p = 0, ±\/3q. 
In Fig. 0(b) we superpose the pattern generated by the 
interference term of the model Wigner function on top of 
the numerically calculated quasidistribution; the agree- 
ment between the two is excellent. Thus, our model ef- 
fectively describes the form of the interference pattern 
and predicts negative areas of the Wigner function. 

Let us emphasize that this analytical model is based 
exclusively on two considerations: the position of the in- 
terfering components in the phase space, and the phase 
relations between them, which were derived from our 
study of the triplet photon statistics for this problem. 
This shows that the interference pattern is very "stiff" , 
i.e. these two considerations strictly impose its specific 
form. Consequently, the interference pattern does not 
change substantially as long as the crucial features of the 



tpi{q) = V2A(-2q) exp(V3ig 2 / 2 - «r/6), 

<p 2 (q) = v / 2-4(-2q)exp(-v / 3^ 2 /2 + i7r/6)- (13) 

Given this result, we can use the approximate form of 
the Wigner function in Eq. (j^) to model the numerically 
calculated Wigner function. Some problems arise from 
the fact that the three components are localized along 
straight lines. In this case the stationary phase approxi- 
mation fails to work for points belonging to the same arm, 
and the Wigner function of each component depends sub- 
stantially on the envelope. Therefore we will denote them 
as W Vo (q, p), W Vl (q,p), and W V2 (q, p) without specifying 
their detailed form. Nevertheless, the stationary phase 
approximation can be safely used to calculate the inter- 
ference pattern between the arms, where the contributing 
points belong to two distinct arms. Thus we represent the 
model Wigner function as a sum of four components 

W v (q,p) = W Va (q,p) + W Vl {q, p) 

+W V2 (q,p) + W iDt (q,p) 7 (14) 

where the interference term Wi n t(q,p) is given by 



(15) 



state remain fixed. In particular, the interaction time 
and the pump amplitude have only a slight influence on 
the basic form of the interference pattern, as they deter- 
mine only the amount of probability density transferred 
to the arms of the quasidistribution. 



V. CONSEQUENCES OF PHASE SPACE 
INTERFERENCE AND EXPERIMENTAL 
PROSPECTS 

We will now briefly review the consequences of these 
phase-space interference effects and the prospects for ex- 
perimental demonstration of quantum interference us- 
ing three photon down conversion. First, let us discuss 
signatures of quantum coherence that can be directly 
observed in the experimental data. An experimentally 
established technique for measuring the Wigner func- 
tion of a light mode is optical homodyne tomography 



A _ 2g _J^U ^L-2q cos M^-Va 



VsJ~\V3 q ' 12 



^ 2 7T 



\p\ < -q 

p > max{ — \/3q, 0} 
p < min{\/3(Z, 0} 
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In this method, the Wigner function is recon- 
structed from distributions of the quadrature operator 
x$ = {ae~ %e + a)e l6 )/^/2, measured with the help of a 
balanced homodyne detector. These distributions are 
projections of the Wigner function on the axis defined 
by the equation qcosO — psin^ = 0. In Fig. || we plot 
the quadrature distributions for the phase 9 in the range 
(0,7r/6). Due to the symmetry of the Wigner function, 
other distributions have the same form, up to the trans- 
formation x — > —x. 

The fringes appearing for x < in Fig. [2] are a clear 
signature of quantum coherence between the two arms 
of the quasidistribution that are projected onto the same 
half-axis. We can describe the position of the fringes us- 
ing the model three-component wave function derived in 
Eq. (|T3|). For simplicity, we will consider only the phase 
6 = 0, for which the fringes have the best visibility due 
to equal contributions from both the arms. The model 
quadrature distribution in the half-axis x < is given by 

\tp!{x) + (p 2 (x)\ 2 = 8A 2 (-2x)cos 2 (V3x 2 /2-tt/6). (16) 

Analysis of this expression reveals some interesting analo- 
gies. Expanding the argument of the cosine function 
around a point x yields that the "local" spacing between 
the consecutive fringes is tt /y/Zx. The same result can 
be obtained by considering a superposition of two coher- 
ent states centered at the points [x, y/3x) and (x, —y/3x), 
i.e., where the contributing pieces of the arms are local- 
ized. Furthermore, the argument of the cosine function 
in Eq. (16) is equal, up to an additive constant, to half 



of the area caught between the two arms of the gener- 
ated state, and the Wigner function of the position eigen- 
state representing the measurement. Thus, the quadra- 
ture distribution given in Eq. ([16]) illustrates Schleich 
and Wheeler's phase space rules for calculating quantum 
transition probabilities |l]]. 

Let us now estimate the effect of dissipation and 
nonunit detector efficiency on the interference pattern 
exhibited by the Wigner function. For this purpose we 
will calculate evolution of the generated state under the 
master equation: 



a) ap — pa^ a), 



(17) 



where 7 is the damping parameter. Evolution over the 
interval At yields the state that is effectively measured 
in a homodyne setup with imperfect detectors charac- 
terized by the quantum efficiency rj = exp(— jAt). In 
phase space, the effect of dissipation is represented by 
coarsening of the Wigner function by convolution with a 
Gaussian function p4[ : 

Wr,{q,p) = ^r I dq'dp'W(q',p') 



2-KT] 



x exp 



1 — T) 

2tj 



This coarsening smears out entirely the very fine details 
of the Wigner function, w hose characteristic length is 
smaller than \f2rj] (1 — 77). In Fig. || we plot the Wigner 
function along the position axis as a function of r;. The 
interference pattern disappears faster in the area more 
distant from the origin of the phase space, where the fre- 
quency of the oscillations is larger. Nevertheless, the first 
negative dip, which is the widest one, can still be noticed 
even for 77 = 0.8. 

Current technology gives some optimistic figures about 
the possibility of detecting the interference pattern, as 
virtually 100% efficient photodetectors are available in 
the range of light intensities measured (in a different con- 
text, that of squeezed light) in a homodyne scheme p5[ . 
However, there are also other mechanisms of losses, such 
as absorption during nonlinear interaction and nonunit 
overlap of the homodyned modes, whose importance can- 
not be estimated without reference to a specific experi- 
mental setup. An analysis of these would be out of place 
here. 

Let us finally consider the impact of pump fluctuations 
on the interference pattern. We illustrate the discussion 
with Fig. f|, depicting the state generated using a noisy 
pump field modelled by a Gaussian P-representation 



P\P) = — exp ' 



(19) 



(18) 



where /3o is the average field amplitude and n is the num- 
ber of thermal photons. In discussing the effect of noise, 
we have to distinguish between phase and amplitude fluc- 
tuations. Phase fluctuations have a quite deleterious ef- 
fect, as a change in the pump phase by $ is equivalent 
to the rotation of the signal phase space by #/3. Con- 
sequently, phase fluctuations average the signal Wigner 
function over a certain phase range. The fringes are most 
fragile near the arms due to neighboring bulk of positive 
probability. The interference pattern in the areas be- 
tween the arms varies slowly with phase, which makes 
it more robust. These properties are clearly visible in 
Fig. || The effect of amplitude fluctuations is not crucial, 
as the position of the fringes does not depend substan- 
tially on the pump amplitude. 



VI. CONCLUSIONS 

We have demonstrated that degenerate three photon 
down conversion generates nonclassical states of light, 
whose Wigner function exhibits nontrivial interference 
pattern due to coherent superposition of distinct phase 
space components. We have developed an analytical 
description of this pattern, which precisely predicts its 
form. Let us note that the rich phase space picture of 
higher order down conversion contrasts with the two pho- 
ton case, where the only signature of quantum coherence 
is suppression of quadrature dispersion P6| . 
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Discussion of the impact of dissipation and pump fluc- 
tuations on the coherence properties of the generated 
state shows that the interference pattern can partly be 
observed even in the presence of moderate amount of 
noise. An important element of the studied scheme is 
that the signal state is generated using a strong exter- 
nal pump, which enhances the usually weak effect of x 
nonlinearity. This allows us the optimism to expect that 
three photon down conversion is perhaps more feasible 
than schemes based on nonlinear self-interaction of the 
signal field. 

The analytical method developed in this paper to de- 
scribe the phase space interference pattern can be applied 
to other cases, where the quasidistribution is a coherent 
superposition of well localized components, for example 
superpositions of two squeezed states [27|, and squeezed 
coherent states for the SU(1,1) group [28|. 



phases equal to multiples of 7r/3. The direction of mo- 
tion can be read out from Eqs. (|a|), showing that the 
sign of the derivative dfl/dr is negative for the phases in 
the intervals (0,7r/3), (2tv/3, 7r), and (47r/3, 5tt/3), and 
positive in the remaining areas. The resulting picture of 
dynamics is presented in Fig. ^. It is seen that the origin 
of the phase space is a three-fold unstable fixed point, 
with the direction of growth 9 = 0, 27r/3, and 47r/3. 



APPENDIX B: ROTATING THE WAVE 
FUNCTION IN PHASE SPACE 

In this appendix we will calculate the rotation of a 
wave function defined by a slowly varying positive func- 
tion A(q) 



Ml) = A{q) 



(Bl) 
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APPENDIX A: CLASSICAL DYNAMICS OF 
THREE PHOTON DOWN CONVERSION 



around the origin of the phase space. An operator per- 
forming this rotation is U(9) — exp(— i9o)a). Its position 
representation is given by 



(q\U(0)W) 



i q + q qq 
■. exp I % 

,-ue\ tan 9 sm9 



(B2) 



where the square root in the complex plane is defined 
by Vre^ = ^fre 1 ^ 12 for r > and -tt < < tt. The 
wave function rotated by an angle 9 is thus given by the 
integral 



The dynamics of multiphoton down conversion under 
classical and quantum equations of motion has been com- 
pared in detail by Braunstein and McLachlan Jig ]; see 
also Ref. p9|| . Here, for completeness, we briefly discuss 
classical trajectories for three photon down conversion in 
the approximation of a constant pump. As the change in 
the pump phase is equivalent to the rotation of the sig- 
nal phase space, we can assume with no loss of generality 
that the pump amplitude f3 is a real positive number. We 
will now decompose the complex signal field amplitude 
into its modulus u and the phase 9. The classical Hamil- 
tonian in this parametrization reads 



H(u,9) = 2\(3u 3 sin 36* , 
and the resulting equations of motion are 



(Al) 



i 



-US 



x / dg'„4(g')exp 



1 q 2 + q' 2 

2 tan0 



. qq 
i — 

sin( 



(B3) 



The stationary phase point for the exponential factor is 
q' = qj cos 9. We will take the value of A(q') at this point 
and perform the integral. Some care has to be taken in 
choosing the proper branch of the square root function, 
when simplifying the final expression. The easiest way to 
avoid problems is to consider separately four intervals of 
9, between 0, tt/2, it, 3ir/2, and 2tt. The final result is: 



,,i(> 



cos 9 V cos 9 



exp 



-— tan 9 
2 



(B4) 



du 2 

— = 3m cos 39, 
At 

— = —3usm39, 
dr 



(A2) 



where r = 2\(3t is the rescaled time. As the energy of 
the system is conserved, trajectories of the system are 
defined by the equation u 3 cos 36* = const. Thus, tra- 
jectories are of hyperbolic-like shape, with asymptotic 
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FIG. 1. The surface (a) and contour (b) plots of the Wigner function representing the signal mode state generated for /3 — 10 
and t = 0.025/A. The dashed lines in the contour plot separate the positive and negative regions of the interference term of 
the model Wigner function derived in Eq. (^|), for comparison with the shaded plot generated from the numerical analysis of 
the full model. 
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FIG. 2. Quadrature distributions {S(x — &$)) for the state 
plotted in Fig. [j]. 
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FIG. 4. The Wigner function representing the signal field 
state generated using a pump with Gaussian noise, charac- 
terized by the parameters fio = 8 and n = 2. The interaction 
time is t = 0.025/A. 
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FIG. 5. Classical trajectories of the signal mode in the 
approximation of a constant pump. 



FIG. 3. The Wigner function along the position axis q after 
dissipation characterized by the parameter f] = exp(— 7 At). 
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